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An attempt is made in this paper to present a concise derivation of the mode theory 
of VLF propagation. Taking note of the fact that the important modes for long-distance 
propagation are near grazing, suitable approximate forms of the wave functions are intro- 
duced at the outset, rather than at the end, of the analysis. It is thus possible to account 
for the influence of earth curvature in a relatively concise manner. The influence of the 
earth's magnetic field is also discussed. Finally, numerical results for the attenuation 
and the phase velocity of the dominant mode are presented. 



1. Introduction 

Despite the extensive work [1-5] 2 on the wave- 
guide mode theory of VLF propagation, certain 
aspects of the subject are not yet resolved. The 
difficulty appears in finding suitable approximations 
to the spherical wave functions in (he rigorous solu- 
tion. It is the purpose of this paper to discuss 
this problem in some detail. For a background of 
the subject the reader is referred to a recent com- 
prehensive paper on the general subject [6]. Here 
the solutions will be obtained from a somewhat 
different point of view. 

2. The Groundwave Field 

We will start by assuming that the source is a 
vertical electric dipole located on the surface of a 
smooth spherical earth of radius a, conductivity a and 
dielectric constant e. Spherical coordinates (r,0,</>) 
are chosen with the dipole located at r=a and 0=0. 
For harmonic time dependence the radial electric 
field component is written in the form 



Erl 



= a(0sin0) 



1/2 



V e i<aT 



(1) 



apart from a constant factor. 3 In the case of an 
airless earth, in which the ionosphere is neglected, it 
is well known that V may be written in the form 
[7,8] 



V^2{wxY' 2 e-" li S 
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where 
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1 Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 

2 Figures in brackets indicate the literature references on page 46. 

3 Equation (1) is normalized such that Vq approaches 2 for /ca->co and A->0. 



and 

\a-\~ieooJ L a^iecoj 

The coefficients t s are solutions of the equation 

w[ (t)-qw! (0 = 0, (3) 

where Wi(t) is an Airy integral and the prime in- 
dicates a derivative with respect to t. In terms of 
Hankel functions of order one third, 



w t (t) = cxp (-2ri/3) (-tt*/3) 1/2 #$[(2/3) (-0 3 



(4) 



The above formula for T T is usually called the 
residue series representation for the groundwave 
field. It is often used as the basis for field strength 
calculations on a spherical smooth earth on the 
assumption that the ionosphere may be neglected 
or separately accounted for. It is valid when 
ka^>^>l and r—a<^<^a. It is also required that the 
values of \t s \ for the important modes are not large 
compared with unity. Such approximations as these 
are certainly valid for VLF waves. 

For the purposes of the subsequent analysis V is 
written as a contour integral in the manner 

V ^e ivl \xl^S *~y *~m dt. (5) 

The contour encloses the poles at t=t s in a clockwise 
sense. The equivalence of eq (2) and (5) is easily 
verified on noting that 



since 



and 



w['(t s ) — (iw' 1 (t s ) = (t s —q 2 )w l (t s ) 
w[(t s ) = qw l (t s ). 
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3. The Sky Waves 

We will now enclose the earth by a concentric 
reflecting shell located at r=a+h as indicated in 
figure 1. The electrical properties of this layer will 
not be specified at the moment, but it is assumed 
that an upgoing wave will be converted to a down- 
going wave. Thus after one reflection 



is converted to 



e ixt Wi{t—y) 



A(t)e- ixt w 2 (t-y), 



where 

w 2 (t) =exp(2^/3) ( -tt//3) 1 /W 1 ^[(2/3) (-*) 3/2 L (6) 

and A{t) is unknown function of t. The boundary 
condition at r=a J rh is written in the form 



c 



'dw(t—y) 



where 



dy *»('-*>] 



y=v 



(7) 



2/o= 



(ka\ 1/s 
— J A*. The quantity A* involves the 

properties of the layer beyond r=a J rh; for the 
moment it is not given explicitly. Formally Ai=Z/rj 
where Z is the radial surface impedance at r=a-\-h. 
Most generally A* (or Z) is a function of t but usually 
it may be taken as a constant [6]. On identifying 
w(f—y) with w 1 (t—y) J rA(t)w 2 (t—y), it readily fol- 
lows from eq (7) that 



U»2(t—yo) + qiW 2 (t—y )J 



(8) 



The downcoming wave characterized by the 
function w 2 (t—y) is now reflected at the ground and 
thus it generates a new upgoing wave of the form 

A(t)B(t)wi{t—y)e- ixt 

The boundary condition at r=a may be written 



B=^«-4--°- 



(9) 



Then on identifying w(t—y) with the sum of the 
downcoming and the upgoing wave it follows that 



B(t)= r^(o-g^2(0 "] 



[(t)-m(tu 



(10) 



The process may be repeated any number of times. 
The resultant field can thus be written as a sum 




<^ggL"r?> 



Figure 1. A sketch of the model employed in this paper. 

The surface of the ground and the lower edge of the reflecting layer are con- 
centric spherical surfaces. 



E,= 



-ikad 



where 



and 



a(0sin0) 



.1=0,1,2 . . . 
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(12) 
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for j even, while 



ixt w 2 (t—y) T 



(t) — qwi(i) 



[A(t)] 2 [B(t)] 2 dt 
(14) 



for j odd. Formally these are geometrical progres- 
sions so they can be summed. On interchanging the 
order of integration and summation, this leads to the 
integral representation 

v rt ,u(x\'* ( 6 e- i * t [w 1 it-y)+A(t)w t (t-y)] /f . 
W J [w' 1 (t)-(iw l (t)}[l-A(t)B(t)] m 



-A(t)B(t)} 

e- i "[w 1 (t-y)+A(t)w 2 (t-y)} 
wi(.t) — quh(t)+A(t)lvi(t) — gw a (t)] 

(15) 



dt. 



Now the contour is to enclose the complex poles which 
occur at t = t n where t n is a solution of 



l-A(t)B(t)=0. 



(16) 



The residue series representation for the total field 
is thus given by 
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V=-2(irx) m e- iT/i 

x y^ e~ i Mwi(t n —y)+A(t n )w 2 (t n —y)] .^ 

As in the residue series formula for the groundwave 
field, eq (17) is valid for fe>>l and h<^<Ca. 

The interpretation of the preceding results is now 
discussed and certain simplifications made to facili- 
tate computation. If (—t)»l, it follows that 

w Ui(o-fl^(oJ 

-L i(-0 1/2 -<z J l ; 

On making the substitution (— £) 1/2 =( — ) 

/ka\ 1/3 
noting that z^=f -— - 1 A, it follows that 



C and 



C+A 



5(f) & 



(19) 



If C is identified as the cosine of an angle of incidence, 
(C — A)/(C+A) can be recognized as a Fresno! 
reflection coefficient. Similarly, if (y — <)>>! 



C — A 



"C'+Ai 



4(f)<co' 



(20) 



who 



/foY' 



(7' = (t/ -0 1/2 . C may be identified 



as the cosine of the angle of incidence at the iono- 
sphere. It is also noted that 



C 



2h\ m 



'-(*+£) 



(21) 



The factor (C'—A,)/(C'+A ( ) is, of course, a Fresnel 
reflection coefficient referred to the bottom of the 
layer at r=a-\-h. 

The modal eq (16) which determines the coeffi- 
cients t n may thus be written in the approximate form 



where 



and 



R g Ri exp (— il)=l=e 
R 



(22) 



<7-A , „ C'-A, 
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= 2 -y[( C2 +W 2 - C3 ] 



h 



if - <C 2 , then S^Stt 1 and £s2fr«7 which cor- 

responds to the modal equation for the flat-earth 
case [5]. 



4. Evaluation of the Residues 

For numerical work it is desirable to express the 
function 



explicitly in terms of Airy functions. 
First, it is noted that 

p V)j(t)-gw 2 (t) l -2i(t n -q 2 ) ( , 

\_dtw[{t)-qw l {t)] t ^ n [w'^tj-qw^tn)] 2 ' K 

where use has been made of the relations 

(1) w"(t) = tw l (t) and 

(2) w 1 (t)w 2 (t)—w[(t)w 2 (t) = —2i which are valid for 
any value of t. Similarly, 



[ 



d w[(t—yo) + q i w 1 (t-y ) l 
dt w , 2 (t—y ) + qiW 2 (t—y )Jt=t n 



2i(t n —y -tfj) 



[^2(tn-yo)+ai^2(t n -yo} 2 



(L>4) 



Thus 



[s>>H,.„ 



-2i{t n -<f) 



[w[(t„)— gwi(t n )][w2(tn)— gw t (t n )] 

2i(t n —yo—qi)[w' i (t n ) — qiV2(t„)] 



(25) 



w'litn) — gwi(t n )][w' 2 {t n — y n ) + qiW 2 (t n — y )} 2 
Therefore, the complete residue series may be written 

V=(irx) m e- ri * /i 



X2 



n = 0, 1, 2| tn~ 



e- ixt n[w 1 (t n -y)+A(t n )w 2 (t n -y)] 

-_£ (t n —y Q —(g)[w 2 (t n ) — qw 2 {t n )} \ 

qw 2 (t n ) [w 2 (t n —yo) + qiW 2 (t n —yo)] 2 J 

(26) 



Numerical evaluation of the series is straightforward 
if tables of the Airy functions of complex argument 
are used. In the situation where 



©v., 



(-«») 1/2 =l 



the Airy functions may be replaced by the first term 
of their asymptotic expansions. Under this approxi- 
mation and after considerable algebraic reduction, it 
is possible to write, for y=0, 



V=±e- i « l \2irkadyi 2 S F- 

w = 0,l,2 . , 



e ikaeC 2 j2 F n 



bt 



% bC 



dC ' BJRi 



c=c n 
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where 









F = 



L„ 



Cl-A* 



and 



^ = 2fo [((7 , )3 _ 
o 



■c* 



(27) 



If the height of the source dipole is z =r Q — a and the 
observer is z—r—a, the F n in the above expression 
should be replaced by 



where 



2j n {z)=B g ~^ exp [ik£ (Cl+^j'dz'j 

+R g ^ 2 exv[~ikf*(0l+^y /2 dz~] 
=E g ^ 2 exp j^[(a 2 +^J /2 ~^]J 

+R^ exp g[(«+f J /2 -^]}- (28) 

If 2zja<^<^Cn it is permissible to use 

2f n (z)^R g ~ 1/2 exp (ikC n z)+R g 1/2 exp (-ikC n z) 

and if R g ^ 1 corresponding to a perfectly conducting 
ground, 

f n (z)Q* cos K7 B 2. 

The value of C n to employ in the residue series is 
related, of course, to t n by 



Hr r 



C n =(^) (-O 172 



K kaJ 



Under the condition that (— if„) 1/2 >>l, the root- 
determining equation becomes 



R t R e ~ i( < 2kJl + v ^ c n~ e -2Tin 

where V is denned by 



(2kh+V)C n =?y i(C' n y-c n *]=t n . 



(29) 



(30) 



Noting that 



Dl 

DC, 



i= 2ka K cn >+W- c °] (31) 

^2khii 2h/a<^C 2 n 
^2k{2ah) ll2 ii 2h/a'>Cl, 



it is seen that the second term within the square 
bracket on the right-hand side is small if kh is reason- 
ably large and if the losses are small. 



5. Discussion of Various Approximations 

As a check on the algebra it may be noted that 
under the flat-earth approximation (i.e., 2h/a<^<^C n 2 ), 
the residue series formula may be written 



T7 2(p/X) 1/2 i ( 2 *p_ ir \ y^ e -i2x P s n /\ j. / n , / n 

y== (A/XT ^ */- A • • • ^ fn(*)M*o), 

(32) 



where 



5n= [- 1+ . atiW/^r 1 



2khR f R g Jc=cn 



P = ad, and S r w =(l-(7 W 2 ) 1/2 ^ 1-^/2. For practical 
purposes, <5 ^l/2 and 8 n ^l for n=l, 2, 3 . . . 
Equation (32) agrees with an earlier analysis [5,6], 
For a first approximation it is probably permissible 
to neglect the curvature altogether. To improve on 
this it is important to retain curvature dependent 
terms in a consistent manner. For example, at VLF, 
the restriction C 2 l ^>2h/a is not as stringent as 

C n >{2lkaY'\ 

Thus eq (27) for V is probably not any better (if as 
good) as eq (32) . To effect a real improvement one 
should retain the Airy functions of argument t n 
without approximation, but allow the Airy functions 
of argument t n —y to be replaced by the first terms 
of their asymptotic expansions. This is justified 
since y is reasonably large compared to unity in the 
VLF range. With this approach the root-deter- 
mining equation may be written 



where 



F.=- 



R i R g F g e- ii2kh+v)c =e- 2 ' rin =l 

w' 2 (t) — qw 2 (t) 
w[{t) — qwi(t) 



[4<-<>»l^£?] 



(33) 



(34) 



i exp | i-H) 



If the ground is assumed to be perfectly conducting 
R g =l and 2=0. Then, on writing 

the mode equation reduces to 



2im=T-2ipC + [(C , y-C i ) =:ka+i arg F g , (35) 



where now 



For long-distance propagation in the VLF range, 
only the mode corresponding to n=l is of importance. 
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Furthermore, since C is small, for a first approxi- 
mation 

w' 2 (t) 



w[(ty 



p iT/G. 



also, is essentially a constant and is approximately 
equal to 1/A*. Thus 



2h\ l/2 



C 2 ^ 



7w 2„ x/2AV /2 _ L0 . i ,/2A\ 



HWM0' 



(37) 



This formula should be a fairly good approximation 
if \C 2 \ is somewhat less than 2h/a. For example, if 
0=0, corresponding to R i = — l ) and if k/a=10~ 2 
the above may be written 



C 2 



0.041248 



-0.01334. 



For h/\ in the region from 2.5 to 5, \C 2 \ is never 
greater than 5X10 -3 . 

Using perturbation methods, higher order approxi- 
mations for the solution of eq (35) may be obtained. 
Details of such computations will not be described 
here. 

6. Certain Extensions of the Theory 

The preceding analysis has referred specifically to 
a vertical electric dipole. If the source is a vertical 
magnetic dipole the results derived above are still 
applicable if the reflection coefficients R g and //, are 
replaced by their counterparts R g and 2?} for hori- 
zontal polarization. The only other change is that q 
is now defined by 

^<7=( — ) — where ~= \ —4 1 • 

Since [q] ^>1, under all practical conditions, 

which is to be used instead of F g . The modes ex- 
cited by a magnetic dipole are of a TE (transverse 
electric) type whereas for an electric dipole they are 
of a TM (transverse magnetic) type. 

7. Application to Specific Ionospheric 
Models 

In the case of a homogeneous isotropic ionosphere 
with a sharp lower boundary, the Fresnel reflection 
coefficient for small C may be written 



n • = r>J — P 

1 C'+A = 



2/SC 



(39) 



where P^-r— In terms of the plasma frequency co 
and the collision frequency v 



]8 ^ (1 — ia)r/u) ( — ica r /o)) ' 



1/2 



where oo r = o) 2 lv under the condition that co<^>. 

Unfortunately, the presence of the earth's mag- 
netic fields renders the ionosphere anisotropic. The 
reflection coefficient is thus a matrix of the form [1] 



[Ri}= 



■R\\ \\R± 
_i_i?n ±R±_ 



where the coefficients ||i?n, \\R± t ±R\\, and ±R± indi- 
cate the complex ratio of a specified electric field in 
the wave after reflection to a specified electric field 
in the wave before reflection. The first subscript 
denotes whether the electric field in the incident 
wave is parallel (n) or perpendicular (j_) to the plane 
of incidence and the second subscript refers in the 
same way to the reflected wave. It is understood, 
of course, that [7i\] is a function of 6 7/ , the cosine of 
the angle of incidence at the ionosphere. Following 
an earner analysis [6], the generalization of eq (33), 
to include the influence of anisotropy, may be 
expressed in the form 

[l-llfiilB.^- 10 * 1 - 100 ] 

X[l-±R±R h g F>~ i(2kh+v)c ] 
-\\R±±R\\R g F 8 R h F h e- i2i2kh + v)c =Q, (40) 

where R g , R h g , F g , F h g , and V have been defined. 
Strictly speaking, equation (40) is only valid when 
there is perfect azimuthal symmetry about the source 
dipole. 

It may be noted that if the conversion coefficients 
IjBjl and _l/?h vanish, the modal equation splits into 
two parts which are 



1 — i\R\\R g F g e 



-i(2kh+V)C_ 



and 



l- ± R ± R*F h ff e- iim+v)c =0. 



(41) 



(42) 



These two uncoupled equations give the solutions for 
the TM and TE modes, respectively. In the more 
general case the modes are coupled. However, for 
low-order modes and particularly for the dominant 
mode (n=l), the coupling is extremely small and a 
good first-order approximation is obtained on the 
assumption that coupling may be neglected. Fur- 
thermore, the TE mode is only of academic or 
secondary interest since it is not excited to any ex- 
tent by a vertical antenna. Thus the constants of 
the first-order mode may be obtained from eq (37) 
with an appropriate value of /3. 
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For propagation in the magnetic meridian or for 
propagation in polar regions the QL (quasi-longitu- 
dinal) approximation [9] to describe conditions within 
the ionosphere is valid. The refractive index jjl, in 
this case, is given by [1] 



fM 2 =l — i(oo r /u) exp (±i<j> L ), 



(43) 



where 



tan0 L 



_o) L longitudinal component of gyrofrequency 

v collision frequency 



and u T = u>l (v 2J rw L 2 ) 1/2 . Thus for highly oblique 
incidence (i.e., C is small) it may be shown that [6] 



where 



P= 



MoMe C +C e 



(44) 



Mo+Me O C e 
where /x and fx e denote the two values of /z, and 

<7o=(1-1/mo 2 )* and C,= (l-l/p.*)l. 

The QL approximation mentioned above is only 
valid when [10] 



4^£ 



<<|(l-1--) 2 

| \ CO CO / 



where u L and u T are the longitudinal and transverse 
components of the (angular) gyrofrequency. As 
mentioned above, this condition is satisfied for 
propagation in the magnetic meridian and/or for a 
steeply dipping field. Clearly, this condition is 
violated when the transverse component of the 
earth's magnetic field is large such as for propagation 
around the magnetic equator. The case of a purely 
horizontal and transverse field has been considered 
by Barber and Crombie [11]. Using their results, 
which are applicable to a sharply bounded ionosphere, 
it may be readily shown that [6] 



\\R\\ 
where now 

and 



o-W 



imR±=0 



(l+i£) 2 -y 2 



(l+iLy /2 (iL-L 2 -y 2 y /2 -iy 



(45) 



(46) 



T CO C0J> , COyCO 

L=— =-=■ and y=— =-• 

co r cog cog 



For east-to-west propagation, y is positive, while 
for west-to-east propagation, y is negative. In the 
present situation it is convenient to define an angle 



<j) T such that 

7=+Z|tan <$> T \ for east-to-west propagation 

= — L|tan <j) T \ for west-to-east propagation. 

For propagation at any angle with respect to the 
earth's magnetic field and for an arbitrary angle of 
dip /, it is necessary to seek more general solutions. 
First, one must solve a quartic equation for the 
refractive indices (ordinary and extraordinary com- 
ponents). Then on matching tangential field com- 
ponents at the air-ionosphere boundary, an explicit 
expression for the reflection coefficient may be ob- 
tained [12]. Using this quite general approach 
Johler and Walters [13] have obtained extensive 
numerical results for the coefficients \\R\\, \\R±, ±R±, 
and j_J?n for various (real) angles of incidence and 
range of values for the electron density, N, collision 
frequency v, dip-angle I and the azimuthal direction 
of the wave <f> a (measured clockwise from the north). 
Formally, the expressions for these reflection and 
conversion coefficients may be substituted in eq (13) 
which, in principle, can be solved for the modes in 
the general case. For VLF and at highly oblique 
angles of incidence (e.g., 82°), the quantity | ni?j_ X 
_l#h| is of the order of 10~ 3 or 10~ 4 and thus coupling 
between the TE and TM modes is negligible. At 
least this is true for the range of values of the iono- 
spheric parameters considered here. Thus the un- 
coupled modal equation may again be used for an 
adequate first approximation after having determined 
the appropriate value of fi. 

8. Presentation of Results 

Using the above approaches numerical results 
have been obtained for the complex values of C for 
the dominant mode. 4 From this quantity, it is a 
simple matter to compute the corresponding value of 
S which describes the propagation in the horizontal 
direction. For example, the attenuation per unit 
distance is proportional to the quantity — H Im S 
and the phase velocity relative to c is 1/Re S. In 
the present paper, these are the quantities of prime 
interest since they are adequate to describe the 
propagation characteristics of VLF radio waves at 
large distances from the source (i.e., >2,000 km). 
At shorter distances higher order modes should be 
considered. Their influence has been discussed 
elsewhere [5]. 

In figure 2 the attenuation coefficient, for the 
dominant mode (n=l), is plotted as a function of 
H(=h/\) for a curved earth where h/a= 0.01. The ( 
upper- and right-hand scales in this figure and in 
many which follow, are a frequency scale and at- 
tenuation scale in decibels per 1,000 km of path 
length applicable to a height of 70 km. The quasi- 
longitudinal approximation is employed and <t> L is 
either 60° or 0°; the latter, of course, corresponds to 
no magnetic field. The constant B on the curves is 



4 For sake of brevity, the subscript 1, to indicate the first mode, is dropped in 
what follows. 
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defined by 
or alternately 



wo 



D L , 1 03 r 

-d =rf where T =— = , 

il 1> co to(z> z +co£ J 17 " 1 



(47) 



5= 



2ttC 



03 r h 

FREQUENCY, kc/s 
10 12 14 16 18 20 22 24 




Figure 2. Attenuation coefficient, — H Im S, as a function of 
H ( = h/\) /or £/ie dominant mode. 
The upper and the right-hand scales are for h =70 km. 

The four B values on the curves (0.02, 0.05, 0.1, and 
0.2) correspond to 03 r values of 13, 5.4, 2.7, and 
1.3 X10 5 , respectively. For a given value of J5, it is 
seen that the attentuation depends on </> L . For large 
values of B (i.e., poorly conducting ionosphere) the 
influence of the terrestrial magnetic field is to lower 
the attenuation whereas the converse is true at 
smaller values of B. Since B is in the range 0.05 to 
0.01 typically for daytime conditions, the influence 
of the magnetic field is small under the quasi- 
longitudinal approximation. 

The relative phase velocity of the dominant mode 
is shown in figures 3a and 3b for <j> L =0 and 60°, 
respectively. The (upper) frequency scale is for a 
height of 70 km. From these it appears that the 
influence of the magnetic field is quite small. It is 
rather interesting to note that for a frequency of 
14 kc/s the phase velocity is very near, c, the velocity 
of light. 5 It should be remembered, of course, that 
the phase velocity is referred to the surface of the 
ground in this particular analysis. 

Since one is interested in relatively small deviations 
of 1/Re (S) from unity, considerable care is needed in 
approximating the wavef auctions involved. The 
results shown in figures 3a and 3b for a curved earth 

5 Recent experimental data appear to confirm this point (M. L. Tibbals, 
private communication). 



are not expected to be accurate to more than a few 
parts in 10 4 . 

FREQUENCY, kc/s 



12 14 16 18 20 22 24 




0.996 



Figures 3a and b. The relative phase velocity, 1/Re S as[a 
function of H for the dominant mode. 

The upper-frequency scale is for h = 70 km. 
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The calculation of the mode characteristics for a 
magnetic field of arbitrary direction is, of course, 
quite involved. However, approximate results can 
be obtained if coupling between TE and TM type 
modes is neglected [6]. Actually, in the case of a 
purely transverse field, the coupling is zero anyway. 

In the following the influence of an arbitary mag- 
netic field is demonstrated by calculating the ratio, 
P, of the attenuation with a magnetic field to the 
attenuation without the magnetic field. Various 
values of the dip-angle /, and the azimuthal 
direction of propagation, <j> a , are chosen. The ratio, 
Q, of the phase-velocity deviations with and without 
a magnetic field is also considered. Specifically, these 
ratios are 



P= 



Attenuation with magnetic field 



Im 8 



Attenuation without magnetic field [Im S] H=0 

(48) 



and 



Q, 



Phase velocity deviation with magnetic field 
Phase velocity deviation without magnetic field 

(49) 

1 

2a 
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To obtain numerical values of P and Q, computed 
data from J. E. Johler [14] were requested for the 
reflection coefficient \\R\\ for a real angle of incidence 
of 82° and at frequencies of 10, 12, 14, 16, 20, and 
22 kc/s. The appropriate value of was then 
calculated from 



\\Ri\=-e-v c ' 



1(50) 



where (7' = cos 82°. The electron density values for 
the sharply bounded model were 10 3 and 3X10 3 
electrons/cm 3 , respectively. The collision frequency 
was taken to be 2X10 7 and the strength of the 
earth's magnetic field was 0.5 gauss. The angle 
<j) a which defined the direction of propagation 
measured clockwise from north assumed the values 
0°, 60°, 120°, 180°, 240°, and 300°. As a check on 
the method, values of /? were calculated directly 
from eq (40) which is applicable to a =9O° and 180° 
for 7=0. 

The results for the attenuation ratios are shown 
in figures 4a and 4b for a dip-angle I of 0°, in figures 
5a and 5b for an I of 45°, and in figures 6a and 6b for 
an /of 84.3°. 

These results would indicate that nonreciprocity 
is most evident in propagation around the geomag- 
netic equator. North-to-south or south-to-north 
propagation appears to be almost reciprocal. It is 
also evident that for this range of electron densities, 
the nonreciprocity for attenuation is more marked 
for the smaller densities. It is also apparent from 
these results that the nonreciprocal effects are 
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Figures 4a and b. The ratio of the attenuation rates with 
and without a (transverse) terrestrial magnetic field. 
4> a is the direction of propagation with respect to geomagnetic north. 
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Figures 5a and b. The ratio of the attenuation rates with and 
without a terrestrial magnetic field for a dip angle of 45° . 
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diminished as the magnetic field approaches the 
vertical. Qualitatively the curves are very similar 
and there is a strong suggestion that nonreciprocity 
is proportional to the component of the earth's 
magnetic field which is horizontal and transverse 
to the direction of propagation. 
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Figures 6a and b. The ratio of the attenuation rates with and 
without a terrestrial magnetic field for a dip angle of 84.3°. 
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Figures 7a and b. The ratio of the phase-velocity deviations 
with and without a (transverse) terrestrial magnetic field. 
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Figures 8a and b. The ratio of the phase-velocity deviations 
with and without a terrestrial magnetic field for a dip angle 
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Figures 9a and b. The ratio of the phase-velocity deviations 
with and without a terresbial magnetic field for a div anale 

of 84.3°. 



The corresponding curves for the ratio Q of the 
phase velocity deviations are shown in figures 7a and 
7b for 7=0, in figures 8a and 8b for 7=45°, and in 
figures 9a and 9b for 7=84.3 °. It is rather interesting 
to see that these are not symmetrical about the 
north-south geomagnetic axis. It is also evident 
that the nonreciprocity in the phase velocit}? 
becomes greater at the lower frequencies. It is 
noticed that Q is never more than 4-percent greater 
than unity. This corresponds to a change of phase 
velocity of only a few parts in 10 5 . 
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9. Concluding Remarks 

In this paper no direct reference has been made to 
experimental data. In another paper [15] the 
experimental data on attenuation rates at VLF in 
daytime [16] were compared with calculations using 
formulas from the present paper. The agreement 
between theory and experiment was quite good and 
the nature of the nonreciprocity between east-west 
and west-east propagation was satisfactorily ac- 
counted for. Unfortunately, the experimental data 
on phase velocity is quite sparse. It would be 
interesting to compare the computed results given 
in this paper with absolute measurements of phase 
velocities at various frequencies over a given path. 
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10.2. List of Symbols 

(r,0,0)= spherical coordinates of observer, 
a = radius of the earth, 
h = height of reflecting layer, 
k = 2 x/ wavelength, 
is central angle, 
x=(ka[2)W 6, 
y=(2/ka)^ (r-a), 
iq=(ka/2)W A, 

A== / iepco y/g ry 16QC0 y/ 2 

W\(t) is Airy function of first type defined by eq (4), 
w 2 (t) is Airy function of second type defined by eq 
(6), 
t is argument of Airy function, 

*7o^l207r, 

Z is surface impedance of ionospheric reflecting 

layer at r = a-\-h, 
%=(k/2) 1/3 A t -, 
y =(2/kay/*kh, 
w(t) is some linear combination of w\(t) and w 2 (t), 
A (t) is defined by eq (8) ; it is related to the iono- 
spheric reflection coefficient, 
B(t) is defined by eq (18) ; it is related to the ground- 
reflection coefficient, 
C=(2Jkay^(-t)V2; Cmay be geometrically in- 
terpreted as the cosine of the angle of in- 
cidence at the ground, 
C'=(C 2J r2h/ay/ 2 ; C" may be geometrically inter- 
preted as the cosine of the angle of inci- 
dence at the ionosphere, 
R g and R{ are reflection coefficients at the ground and the 
ionosphere, respectively, 

V is a dimensionless quantity proportional to the 

vertical electric field, 
Vj is defined by eq (13) ; it may be interpreted 

loosely as the j'th order sky wave, 
t n is the n'th root of eq (16); it determines the 
attenuation and phase of the n'th mode, 
p = ad; it is the great circle distance between 

source and observer, 
8 n is defined by eq (32) ; it is a coefficient in the 
mode sum, 
F g is defined by eq (34) ; it is a wave correction to 
the ground-reflection coefficient, 

W r = 0)o 2 Jj>, 

coo is the (angular) plasma frequency, 
v is the collision frequency, 

j8 is defined by eq (39) ; it characterizes the iono- 
spheric reflection coefficient \\R\\, 

[Ri} = \ n n ; this is a matrix which describes 
L_i_it|| ±Rlj 

reflection from an anisotropic layer, 
/x is defined by eq (43) ; it is a double- valued index 
of refraction for the ionosphere subject to 
the QL approximation, 
juo and fi e are the two components of /*, 

co L is the longitudinal component of the gyrofre- 

quency, 
cot is the transverse component of the gyrofre- 
quency, 

B=L/H=^> 

co r /i 
T co co^-fcoz, 2 ) 1 ' 2 



co r co 

0L=tan- 1 (co L W, 
P and Q are defined by eqs (48) and (49), 

7= dip angle of the earth's magnetic field, 
cf> a is the azimuthal direction of propagation of the 
wave; c> o =0 corresponds to propagation 
from south-to-north, <£ a = 90° corresponds 
to propagation from west-to-east, etc. 
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